When Do Internal Shocks End and External Shocks Begin? 
Early-Time Broadband Modelling of GRB 051111 



N. R. Butleri'2, W. Li^, D. Perley^, K. Y. Huang^ Y. Urata6'^ J. X. Prochaska^, J. S. 
Bloom^, A. V. Filippenko^, R. J. Foley^, D. Kocevski^, H.-W. Chen^ Y. Qiu^ P. H. Kuo^ 
F. Y. Huang^ W. H. Ip^ T. Tamagawa^ K. Onda^ M. Tashiro^ K. Makishima*^'^ S. 
^ ' Nishihara^", and Y. Sarugaku^° 

o ■ 

O 
(N 



> 
O 

m 

m 
> 



o 
o 



X 



ABSTRACT 

Even with the renaissance in gamma-ray burst (GRB) research fostered by 
the Swift satelhte, few bursts have both contemporaneous observations at long 
wavelengths and exquisite observations at later times across the electromagnetic 



, spectrum. We present here contemporaneous imaging with the KAIT robotic 

! optical telescope, dense optical sampling with Lulin, and supplemented with 

o . 

\0 . For the first time, we can test the constancy of microphysical parameters in the 



internal-external shock paradigm and carefully trace the flow of energy from the 
GRB to the surrounding medium. KAIT data taken ^ 1 minute after the start 
of GRB 051111 and coinciding with the fading gamma-ray tail of the prompt 
emission indicate a smooth re-injection of energy into the shock. No color change 
! is apparent in observations beginning ~ 1.5 minutes after the GRB and lasting 

for the first hour after the burst. There are achromatic flux modulations about 
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the best-fit model at late {t ^ 10^ s) times, possibly due to variations in the 
external density. We find that the host-galaxy extinction is well fit by a curve 
similar to that of the Small Magellanic Cloud. Low visual extinction, Ay f« 
0.2 mag, combined with high column densities determined from the X-ray and 
optical spectroscopy {Nh > 10^^ cm~^), indicate a low dust-to-metals ratio and 
a possible over- abundance of the light metals. An apparent small ratio of total to 
selective extinction {Ry ^ 2) argues against dust destruction by the GRB. Time 
constancy of both the IR/optical/UV spectral energy distribution and the soft 
X-ray absorption suggests that the absorbing material is not local to the GRB. 

Subject headings: gamma rays: bursts — supernovae: general — X-rays: general 
— telescopes 



Introduction 



The Swift gamma-ray burst (GRB) satellite (IGehrels et al.ll2004l ) continues to unleash a 
torrent of finely time- and energy-sampled photons arising from the bursts and their impact 
on the surrounding medium (i.e., the GRB afterglows). The X-ray Telescope (XRT) provides 
detailed light curves for 1-2 bursts and burst afterglows per week, allowing us to routinely 
view the practically uncharted first 100s to ~ 1 day in the life of a GRB. In the X-ray data, 
we see complex phenomena: unexpected rapid decays, also sometimes temporally flat an d 
apparently re-energized afterglows, and massive X-ray flares (see, e.g., iNousek et al.ll2006l ). 
Unfortunately, it has proven very challenging to complement these data with an early and 
rapid-cadence data set at longer wavelengths. Because studies of GRB afterglows prior to 
Swift rely primarily on observations in the optical and longer wavelengths, and mostly at 
later times, contemporaneous observations are critical for connecting new-found insights to 
the large body of previous work. Correlated, broadband observations of events which exhibit 
extremely energetic broadband emission from the radio to the gamma-rays (e.g.. Figure [3]) 
allow us to study and potentially understand the new phenomenology and to pose answers 
to numerous open questions. 

Optical observation of GRBs in the prompt phase and shortly there after are extremely 
rare. Of ~ 50 optical afterglows detected prior to Swift, only one (GRB 990123. lAkerlof et al. 



19991 ) was detected during the prompt phase and only a handful were de tected in the first 
10 min of the afterglow: GRB 021 004 JFox et all lioOSal) , GRB 021211 jibx et al.lboOSbl : 
Li et al.l[2003bl : IVestrand et al.ll2004h . and GRB 030418 jRvkoff et al.ll2004[ ). In the Swift era, 
the rate of early detections has jumped dramatically, thanks to rapidly communicated and 
tight GRB localizations from Swift and also due to the maturing system of ground-based 
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robotic telescopes (e.g., ~15 early detections by ROTSeI"] alone). The UV Optical Tele- 
scope (UVOT) on Swift has the potential to match this performance, particularly with 
a recent prioritization of early unfiltered observations of the very red afterglows. Perhaps 
most impressi ve, the RAPTOR exp eriment has detected t wo GRBs during the prompt phase: 
GRB 041219 JVestrand et allboosh and GRB 050820A, JVestrand et al.ll2006h . The prompt 
and contemporaneo us long wavelengt h emission of GRB 041219A was actually discovered 
at IR wavelengths (IBlake et al.l |2005| ) with the Peters Automated Infrared Imaging Tele- 



scope ( iBloom et al. 



20061 ). Even though the Galactic extinction toward the GRB was large. 



GRB 041219A is the only prior burst with prompt, long-wavelength observations at multiple 
frequencies. 

Here we present observations of GRB 051111 conducted with the robotic 0.7 6-m Katz- 
man Automatic Imaging Telescope (KAIT; Filippenko et al. 2001 : Li et al. 2003al lbl). begin- 
ning with an unfiltered exposure 43.7 s after the GRB trigger from the Burst Alert Telescope 
(BAT) on Swift and just catching the tail of the prompt emission. We obtained color in- 
formation in the form of V and J-band observations beginning just 73.7 s after the trigger. 
Supplemented with B, V, R, and /-band observations with the Lulin One- meter Telescope 
and J, H, and Kg-hand observations with PAIRITEL taken hours after the burst and with 
other observations reported in the literature, we present an unrivaled and impressive broad- 
band view of this afterglow and its early evolution. 

The study of the early emission during and following a GRB is a critical step toward un- 
derstanding the origin of the emission and the nature of the surrounding medium. It is widely 
acce pted that GRBs are produced by a self- interacting relativistic outflow ( "internal shocks," 
e.g.. lFenimore. Madras, fc Nayakshinlll996l ) which heats and shocks the sur rounding medium 



( "external shocks" ) to for m a long-lasting afterglow at longer wavelengths (IMeszaros fc Rees 
19971 : ISari fc PiranI Il999l ). Transient emission from a "reverse shock," which propagates 
backward toward the central engine in the shock frame, can also be produced. Despite 



an early indication to the contrary from GRB 990123 (e.g., lAkerlof et al.l Il999l ). the re- 
verse shock does not appear to b e a common feature in the early optical data (see, also, 
McMahon. Kumar. &: PiranI |2006| ) . Primarily due to the larger physical size, the external 
shocks are expected to generate smoother light curv es, and the ^ 10% of GRBs with smoot h 



time histories may require only external shocks (e.g.. lMcMahon. Kumar, fc Panaitescull2004j ). 



However, the growing sample of early afterglow observations militates against such a 
simple interpretation. As mentioned a bove, the X-ray data (and also some optical data, e.g. 



Fox et al.l l2003al : IWozniak et al.l |2005| ) evidence shock refreshment (IRees fc Meszarod Il998 



11 
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2000 
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■lamirez-Ruiz. Merloni. & Ree 


3 2001 


) or continued central engine ac- 


tivitv ( 


Rees & Meszaros 


2000 


MacFadven et al. 


2001; 


Ramirez-RuizlboO^ 


:lLee & Ramirez- Ruiz 


2002 


), which can persist for several hundreds of seconds or longer (e.g., 


N^ousek et al. 


2006 


). 



While its physical size is still small compared to the timescale for synchrotron cooling, the 
external sh ock can cool rapidly. Energy can be efficiently radiated from the shock at this 
stage (e.g., ISari. Piran. fc Narayanlll998l ). prior to the onset of an adiabatic evolution. De- 
pending on the breadth of the shock shell, which depends on the duration of the period of 
internal shocks generation, the internal shock emission can overlap with the external shock 
emission. Below, we use the exquisitely sampled data for GRB 051111 to disentangle the 
early-time emission components and also to study the dust and gas properties of the host 
galaxy. 



Observations and Data Reduction 



At 05:59:41.4 UT, Swift BAT triggered and localized GRB 051111 (ISakamoto et al. 



20051 ). The 3' radius error region fr om BAT enabled the ra pid detection of a new optical 
source only 26.9s after the trigger (iRujopakarn et al.l l2005l ). while the burst was still in 
progress. A barrage of detections in the optical, ultraviolet (UV), infrared (IR), radio, and 
X-ray bands followed. Figure 1 shows a snapshot of the full data set, interpolated via the 
modelling below to t = 100 s. 



2.1. Gamma- and X-ray 

We downloaded the Swift BAT gamma-ray data from the Swift Archivj^. The energy 
scale and mask weighting were established by running the bateconvert and batmaskwtevt 
tasks from HEAsoft 6.0.4. Spectra and light curves were extracted with the batbinevt task, 
and response matrices were produced by running batdrmgen. We apply the systematic error 
corrections to the low-energy BAT spectral data as suggested by the BAT Digest websitJ^. 



and fit the data using ISII 



J 14 



The Swift XRT X-ray data from 5 follow-up observations of the GRB 051111 field 
were downloaded from the Swift Archive and reduced by running version 0.10.3 of the 



^^ftp:/ /legacy.gsfc. nasa.gov/swift /data 

^■^http: / / swift.gsfc.nasa.gov/docs / swift /analysis /bat jdigest.html 
i4http://space.mit.edu/CXC/ISIS/ 



- 5 - 



xrtpipeline reduction script from the HEAsoft G.O.q^j software release. From there, we 
bin the data in time, exclude pileup chip regions for each time interval, ac count for los t 



flux due to bad pixels, and produce spectra using custom IDL scripts (e.g., iButlerl |2006| ). 
The data cover the time range from 5.56 ksec until 405.29 ksec after the burst, with a total 
exposure (livetime) of 51.98 ksec. Spectral response files are generated using the xrtmkarf 
task and the latest calibration database (CALDB) files (version 8, 2006-04-27). The spectra 
are modelled using ISIS. For each spectral bin, we require a signal-to-noise ratio (S/N) of 
3.5. We define S/N as the background-subtracted number of counts divided by the square 
root of the sum of the signal counts and the variance in the background. We define the 
background region as that where the number of counts in an aperture the size of the source 
extraction region is within 2a of the median background over the chip in that aperture for 
one contiguous follow-up observation. For the Photon Count (PC) mode data, the source 
aperture is a circle of radius 16 pixels. 



2.2. Optical 



The robot ic 0.76-m Katzman Automatic Imaging Telescope (iFilippenko et al.l 12001 
Li et al.l[2003al lbl. KAIT;) at Lick Observatory observed GRB 051111 in a seri es of images au 
tomatically obtained starting at 05:60:25.2 UT (43.7 s after the BAT trigger: Eiet al.ll2005[ ). 
The sequence includes a combination of images taken with the V and / filters, as well as 



some that are unfiltered (Table 1). The optical transient first identified by lRujopakarn et al. 
J2005h at a = 23'*12'"33.2^ 5 = +18°22'29.1" ('J2 000.0) was clearly detected in each ex- 



posure. We began B, V, R, and / band imaging (jl; 



uang et al.ll2005h of the GRB 051111 



afterglow using the Lulin One-meter Telescope (LOT; iHuang et al.ll2005bl : lUrata et al.ll2005l ) 
at 10:07 UT. We detected the optical afterglow clearly in each band until observations ceased 
at 14:42:06 UT (Table 2). 

We find that the combination of the KAIT optics and the quantum efficiency of the 
Apogee CCD camera makes the KAIT unfiltered observ ati ons mostly mimic the R band. 
We determine a small color correction ( Li et al. 2003a| jbl: iRiess et al.l Il999l ) between the 
unfiltered and R band photometry of 0.11 ± 0.02 mag using the {V — R) and {R — /) color 
information from the Lulin observations. 

We use point spread function (PSF) fitting through the IDL DAOPIIOT0 package to 



^'"'http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
^®http://idlastro. gsfc.nasa.gov 
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reduce our data. Only unsaturated and spherically symmetric sources within a given CCD 
exposure are used to model the stellar PSF for that ex posure. The ab solute photometric 
calibration for KAIT is determined using more than 10 iLandoltl (Il992l ) fiel ds observed on 
Nov. 24 UT at a range of airmasses. The Lulin photometry is calibrated using iLandoltl (Il992l ) 
fields PG0231+051, SA92, and SA97, and also cross-checked using the KAIT standard star 
observations. 



2.3. Near-IR and UV 



After initially poor transmission conditions (i.e., clouds) at Mt. Hopkins, we obtained a 
total imaging exposure of 188 s on the field of GRB 051111 before hitting a telescope limit. 
In sim ultaneous observa t ions w ith PAIRITEL (IBloom et al.l 120061 ) we measure magnitudes 
of the iRuiopakarn et all f 20051) transient of J = 16.55 ± 0.03, H = 15.85 ± 0.04, and = 
15.29 ±0.06 Jeioom et al.ll2005h . relative to 2MASS (Two Micron All Sky Survey). The data 
were taken from 07:23:29 to 07:28:26 UT, centered 1.44 hrs after the BAT trigger. 



We refine the Swift UVOT photometry (jPoole et al.ll2005l ) in th e V, B, and U bands by 
employing a tight extraction region in order to maximize the S/N (see lLi et al.ll2006l ). We find 
V= 19.38±0.15 {t = 16298 -17198 s post BAT trigger), B = 20.04±0.10 {t = 11237-12136 
s post BAT trigger), and U = 19.81 ± 0.12 (t = 10329 - 11229 s post BAT trigger). These 
numbers are consistent with those of iPoole et al.l (120051). but with considerably srnaller error 
bars. We extrapolate the Galactic extinction curve ( Schlegel. Finkbeiner. fc DavisI 19981 ) into 
the UV in order to perform the fitting for the UVWl, UVM2, and UVW2 bands below. 



3. Fits to the Data 

3.1. Optical Light Curve and Lack of Color Change 

As shown in Figure (H the KAIT R band data are well fit by a broken powerlaw with 
tbreak = (700 ± 260) s {x^ / u = 5.61/11). The best-fit flux decay indices are shown in Figure 
[U The fit is statistically unacceptable without the break ix^ /f = 65.15/13, a = 0.83±0.01). 
Including the Lulin R band data in the fit, the model parameters do not change, but the fit 
quality degrades (x^/z^ = 38.9/20). Similar fit qualities are found for the data in the B, V, 
and I bands by applying magnitude offsets to the best-fit R band model (x^/j^ = 27.0/11, 
54.1/13, and 64.0/10, respectively). From the fits, we derive the following afterglow colors: 
{B-R) = 0.97 ± 0. 02, (V-R) = 0.37 ± 0.01, and (i? - /) = 0.80 ± 0.01 mag, corrected for 



Galactic extinction (jSchlegel. Finkbeiner. fc DavisI Il998l ) 
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The [R — I) and {V — R) colors, without the Galactic extinction correction, are plotted 
in the sub-panel of Figured], and the multi-band data are plotted in Figure [21 along with the 
fits. The KAIT and Lulin data are consistent with no color evolution, and the light-curve 
break is most likely achromatic. From the {V — I) colors before and after the break, we 
derive a change in the spectral index 6(3 < 0.38 (3(t). Figure [2] displays another interesting 
characteristic of the data: there are strong residuals in the Lulin data with St/t ~ 0.1 relative 
to the broken powerlaw fit, and these are correlated across spectral bands. 

As we discuss further below, the PAIRITEL data, as well as the available data from 
the GCN in the optical and UV, are consistent with the broken powerlaw decay, modified 
by dust absorption. 



3.2. Gamma- and X-ray 



The BAT gamma-ray light curve exhibits a characteristic FRED-like ( "fast rise exponen- 
tial decline") time decay, with duration Tgo ~ 30 s. From a wavelet analysis, we determine 
the start of the burst as 7.17 s prior to the BAT trigger (Figure [6]). From this point onward, 
we use this time as the beginning of the burst. This has the effect of increasing the initial 
optical decay index from 0.77 to 0.81 (Figure [I]), but does not affect any of the other results 
above. The 15-150 keV spectrum from the beginning of the burst to 40 s later is adequately 
fit by a powerlaw with photon index F = 1.27 ± 0.05 {x^ /v = 53.65/51). The model fluence 
in the 15-1 00 keV band is (2.1 ± 0.1) x 10~^ erg crn"^. T he fits are consistent with those 



reported by lSakamoto et al.l (l2005l ) and lKrimm et al.l ((20051), where a modestly longer burst 
time extraction region is applied. 

There is a significant amount of spectral evolution, comparing the time prior to burst 
peak to the time post peak. For the burst rise from s to 10 s, we find Fi = 1.11 ±0.07 and 
^iskev = (0.48 ±0.04) mJy {x^/u = 49.35/51). For the burst decline from 10 s to 40 s, we find 
T2 = 1.38±0.08 and Figkev = (0.31±0.03) mJy (x'^/u = 41.15/45). The photon index for 20 
s to 40 s (Fa = 1.4±0.1, x^/^^ = 29.28/31) is consistent with the 10 s to 40 s value, indicating 
that this value provides a reasonable characterization for the full dec line phase. The burst 
was also observed by the Suzaku Wideband All-sky Monitor (WAM; Yamaoka et al.l 12^006) 
and found to have a similar time profile as from Swift and also a consistent time integrated 
spectr al slope (F = 1.5 ±0.3 ) and a fluence in the 100-700 keV band of (8.4 ±0.8) x 10"^ erg 



cm ^ (lYamaoka et al.l |2005| ) . Combining the fluence determinations from Swift, and Suzaku, 



we find a 15-700 keV fluence of (1.05 ± 0.08) x 10~^ erg cm~^. Because the peak vFy energy 
for this burst has not been measured, the bolometric fluence could be substantially larger 
than the 15-700 keV fluence. Assuming that all the prompt photons have been accounted 



- 8 - 



for, a lower limi t on the isotropic equivalent energy emitted in the host frame at z = 1.55 
( iHill et al.ll2005l ) is E[so = 6.2 x 10^^ erg. Here and throughout, we assume a cosmology with 



{h, ^a) = (0.71, 0.3, 0.7). 

The X-ray light curve in the 0.3-10.0 keV band is reasonably well fit {x^/ v = 38.7/32) by 
a powerlaw model with a = —1.6 ±0.1. The spectrum in the same band contains 573 counts 
and is well-fit ix^ = 43.61/36) by an absorbed powerlaw with photon index T = 2.3 ±0.2, 
absorbing column = 1.8 ± 0.4 x 10^^ cm~^, and unabsorbed flux (5.0 ± 0.5) x 10^^^ erg 
cm~^ s~^, for the time period 5.56 to 405.29 ksec after the GRB. We use the corresponding 
average spectral flux at 1 keV, (6.6 ± 0.9) x 10~^ /xJy, to translate the XRT count rate to 
/xJy below. The absorption i s greater than the inferred Galactic value in the source direction 
(^//.Galactic = 5 X 10^° cm'^. IPickev fc Lockmarj|l990l ) a t 4.1 cr significan ce (Ax^ = 16.48, for 



1 additional degree of freedom). At redshift z = 1.55 (iHill et al.ll2005l ). the excess column 
is (7 ± 4) X 10^^ cm~^, allowing for a 20% uncertainty in the Galactic column. We find no 
evidence for ;^ la significant er nission lines in the spec trum. These fits are consistent with 
the preliminary fits reported by iLa Parola et al.l (120051 ). 



4. Discussion 



4.1. Late-Time Broadband Afterglow Modelling 



After t ^ 10^ s , the broadband afterglow data are well described by the external shock 

synch rotron model JPaczvhski fc Rhoadslll993l : lKatzll 19941 : IWaxmanlll997l : IWijers. Rees. fc Meszaros 
19981 ). We do not attempt to model the slowly decaying early optical light curve or the achro- 
matic (or nearly so) optical break discussed above, nor do we fit for the BAT data or the 
apparent fluctuations in the optical at t ~ 6000 s (Figure [2]) at this stage. These features 
are discussed below. Due to the fluctuations in the optical about the best flt model at 
t ~ 10^ s, we are also forced to add a 10% systematic uncertainty component in quadra- 
ture to the statistical uncertainty of each data point, in order to obtain an acceptable flt 
(x^/u = 118.1/103, t > 700 s). Figure S shows the best-flt model, plotted over the radio, 
optical, UV, IR, and X-ray data at t > 700 s. In the next section, we describe the modelling 
of the host-frame extinction by dust, which has been taken into account (Figure E]) in the 
figure so that the data in the optical, UV, and IR appear on one common curve. 

Because the optical time decay is shallow and the X-ray decay is steep, we expect and 
observe the data to be well fit by a uniform density, interstellar medium (ISM) model. A 
model with n oc R~^, describing the expected density contours for a progenitor star with a 
significant wind, would exhibit a synchrotron cooling break frequency Uc which increases in 
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frequency with time (IChevalier &: Lil 120001 ) . This leads to shallower decays at high versus 
low energy, opposite to the behavior of the ISM model. For the const ant- density model, the 
optical and X-ray decay indices (ctopt = — 0.96±0.04, ax = — 1.6±0.1) and the X-ray energy 
index {(3x = —1.3 ±0.2) constrain the powerlaw index p describing the energy distrib ution of 



synchrotron-emitting electrons: p = 2.35 ± 0.05 (e.g., iSari. Piran. fc Narayanlll998l ). In the 
best-fit model, the synchrotron cooling break lies between the optical and X-ray bands during 
the observation. The predicted value for the optical energy index is Popt = —0.68 ± 0.03. 



The model contains 4 additional free parameters (e.g., ISari. Piran. &: Narayanl Il998l ): 
7]^, ee, es, and n, the efficiency factor relating the shock energy to the energy released in 
gamma-rays, the fractions of equipartition energy going into the electrons and magnetic field 
B, and the density n. These are constrained by the three observed radio, optical, and X-ray 
fiuxes and by the equation describing the time evolution of the synchrotron spectral peak 
frequency Um- If we assume that the peak in the synchrotron spectrum passes through the 
optical band at the same time the light curve breaks {t ~ 700 s), as might be expected if the 
early optical light curve is dominated by the internal shocks, then we derive the parameters 
listed in Table 3. As we describe below, the pre-break optical light curve may also be due 
to external shocks, in which case we only have an upper limit on the passage time tm.opt for 
the synchrotron peak through the optical bands. Aside from e^, the model parameters are 
sensitive to (and roughly linearly proportional to) tm,opt (Table 3). 



4.2. Optical/IR/UV Spectral Energy Distribution 



A tight constraint can be placed on the optical absorption considering only the observed 
optical/X-ray spectral fiux ratio -Fj,=^ obs/-^i/=i kev = 467 ±2, at t = 2 x lO^s in the observer 
frame. Here, the X-ray fiux is corrected for absorption. Allowing for the possibility that the 
synchrotron cooling break (z/c) is between observer frame R ban d, un^chi^, and 1 keV in t he 



observer frame, ui kev, the observed optical to X-ray fiux is (e.g., iGalama fc Wijersll200ll ) 



Fu=R,ohs/ ^"^=1 keV 



T^l keV 

^R,obs 



l/2-/3opt 



-1/2 



10 



-OAAt 



where A/j^obs is the absorption at R band in the observer frame. The constraint on jSopt 
derived above for the const ant- density model leads to Aji^ohs < 1-35 mag. For the wind 
density medium, there cannot be a spectral break due to the rapid X-ray versus optical light 
curve decline, and we have A^^ohs = — 6.8/5opt — 6.7. This relation can only be satisfied with 
positive absorption if —0r,pt > 1- This is j ust possible in the wind scenario for < I'R^ohs, 
implying p = 1.6±0.1 (IChevalier &: Lill2000l ) and Aji ^hs < 0.1 mag. The rest-frame extinction 
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in the V band (i.e., 5556 A), Ay, will be a factor of a few times smaller than this, depending 
on the extinction law. 

Figure Eti' displays extinction-law fits to the optical/IR/UV spectral energy distribution 
(SED) interpolated to 5200 s. We interpolate all data using the broken powerlaw fit de- 
termined for the unfiltered KAIT and R band Lulin data (Figured]). We choose 5200 s so 
that little interpolation is required for the IR or UV data. Because there is a prominent 
departure in the IR, optical, and UV light curve data from the fit curve in the time interval 
5 X 10^ — 1.2 X lO^s in Figure H] (see, also. Figure [1]) — which can be modelled with a single 
magnitude offset of ~ 0.6 mag from the powerlaw curve — we must include a variable offset 
in the SED fitting to describe the fiux offset in this time interval. (Explanations for this be- 
havior are discussed below.) Also, we use the constraint derived above for the (unabsorbed) 
optical energy index Popt- 

We begin by fitting the empirical Milky Way, Large Ma gellanic Cloud (LMC), and 
Small Magellanic Cloud (SMC) extinction curves of iPeil (119921 ). The fits (Figure [5^) yield 
rest-frame V band extinction values of Ay = 0.38 ± 0.15 {x^/u = 12.23/7), 0.35 ± 0.08 
{x'^/jj = 5.96/7), and 0.23 ± 0.07 {x^/u = 4.32/7) mag, respectively. The excellent fit of 
the SMC model and the poor fit of the Galactic extinction model point toward a weak 
2175 A dust "bump" and an increased level of far-UV extinction relative to that found in 
the Galaxy. Add i tional ly, we fit the data using the general, 8-parameter family of models 
from iFitzpatrickl (119991 ). focusing on variations in Ay and the ratio of total to selective 
extinction, Ry = Ay/E{B — V). In order to fold in prior knowledge of physically realistic 
values for Ry while not suppressing possible variations in the other parameters describing 
the extinct ion curve (e.g., the magnitude of the dust bump, which we marginalize over), we 
exploit the iReichartI (120011 ) prior. The best-fit extinction curve has Ay = 0.2 ±0.1 mag, in 
agreement with the SMC value determined above, and Ry = 2.0 ± 1.0 ijd 1'^ = 4.9/6). The 
model extinction at R band in the observer frame is A/j obs = 0.7 ± 0.1 mag, consistent with 
the value derived a bove from the optica l/X-ray fiux ratio alone. Finally, we note that the 
extinction curve of ICalzetti et al.l (120001 ) provides a poor fit {y^jv = 9.42/6, with best-fit 
Ay = 0.24 mag, Ry = 2.0). 



Stan dard Galactic dust has Ry = 3.1 (ISneden et al.lll978l ). and the low inferred Ry (see 



also, e.g.. lKrisciunas et al.ll2000l . for the appearance of such values in the SEDs of supernovae) 
implies a dust-grain distribution skewed toward fine grain sizes. The posterior probability 
contours in Figure [5]d show Ay < 0.45 mag {2a). Larger values of Ay /3opt ~ 3 mag are 
possible if we relax the constraint on Popt and allow positive values, however, these models 
are inconsistent with the temporal decay observed in the light curve and cannot produce 
the observed optical/X-ray fiux ratio. For the wind-density medium, which is excluded by 
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the rapid X-ray versus slow optical flux decay unless the early X-ray light curve is actually 
dominated by low-level flaring, the spectral slope can be as shallow as jSopt = —0.3 (z/c above 
the X-ray band, in conflict with the constraint derived above from the optical/X-ray flux 
ratio). But this has little effect on the inferred Ay, and we flnd Ay < 0.5 mag {2a). 

Both the SED flts and the optical/X-ray flux ratio demand a considerably smaller dust 
column than we infer from the soft X-ray absorption (Section l3.2l) or fror n detailed modelling 
of th e GRB 051111 optical spectrum. The Galactic — Ay relation (jPredehl fc Schmitt 



19951 ) and the observed X-ray column irnply Ay = 4 ± 2 mag. (One caveat here is that XRT 
low energy calibration efforts are ongoing^.) This indicates that the ISM exhi bits a low dust- 
to-gas ratio, si milar to the SMC which is ~ 8 times smaller than Galactic (iBouchet et al. 



19851 : |Peilll992l ). Another contributing factor could be a n overabundance of the light me tals 
(which dominate the opacity at soft X-ray wavelengths; Morrison &: McCammon Il983h. as 



also suggested by optical spectroscopy in the case of GRB 050401 (IWatson et al. 



Similar to GRB 050401, the Keck/HI RES spectrum of GRB 051111 flProchaska et al. 



2006). 



2005 



Penprase et al.l l2006l : iProchaskal |2006| ) exhibits strong absorption lines in the trace element 
Zn (implying log(A^H,gas) = 21.2 ±0.2, for solar metalicity), which combines with the inferred 
light metal column from the X-ray data to imply a factor of 4.0^2 e overabundance in the 
multiple-a eleme nts relative to the Fe -group elements. Because of the saturation of the Zn 
absorption in the IWatson et al. J2006h GRB 050401 spectrum and because Zn is only a trace 
element and does not contribute directly to the opacity (e.g., at X-ray wavelengths), this 
possibility should be approached cautiously. 

In a separate paper, we plan to study in more detail the X-ray absorption in this burst 
relative to the absorption properties determined from the optical spectroscopy. From the 
metal abundances inferred from the GRB 051111 Keck spectrum (e.g., Si/Fe or Zn/Fe), we 
flnd evidence that the gas is highly depleted by dust. There is also marginal evidence that 
the light metal column (in gas form, \og{N = 21.0^q;5, from log(A^5j+) = 16.6^q;5) may 
be a factor ~ 10 smaller than that inferred from the X-ray data. It is not clear how these 
column densities are to be reconciled with the constraints imposed by the optical and X-ray 
SED and the soft X-ray absorption. The Keck spectrum shows, for example, Nsi++ < Nsi+, 
which argues against an explanation in terms of a signiflcant ionization of the column. The 
large X-ray column is probably not signiflcantly due to intervening systems (a Mgll absorber 
in the optical spectrum at z = 1.189 only contributes log^Nn) ~ 19.5). 
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4.3. External Shock Origin for the Early Optical Light Curve and GRB? 



From the optical colors, we derive above a stringent constraint on the change in the 
spectral slope across the break at t ~ 700 s, which excludes an explanation based on a 
passage through the spectral bands of the synchrotron cooling break (yielding A/3 = 0.5; e.g., 



Sari. Piran. fc Narayanlll998l ). The most commonly invoked cause of achromatic breaks- 



flux decremen ts due to finally viewing the edge of the sideways expanding GRB jet (e.g., 
Rhoadslll999l ) — is ruled out here by the mildness of the break (Aa = 0.19 ± 0.04) and by 



the extremely slow optical fade after the break. Another possible explanation, which grafts 
directly onto the t > 700 s solution discussed above, is that the density prior to the break 



is inc reasing with radius. A decreasing density, as for a wind medium (e.g., iRykoff et al. 



2004 ). does not work because the fiux increase due to the decreasing optical absorption is 



over-compensated by the decrease in synchrotron fiux with radius. We would also expect to 
see a strong color change due to the changing absorbing column. For the increasing density 
model, the absorbing matter close to the burst is unimportant, and there is little expected 
evolution in the absorption with time. 

The pre-break decay is shallower by Aa = s/(8-|-2s), for a density rising as n cx i?* and 
s = 4 ± 2. This is because the fiux between the synchrotron peak frequency and the cooling 
break is proportional to the sy nchrotron peak frequency z/m"^"*^^ times the peak fiux F„ 



Sari. 



'iran. fc Narayan 



le.g. 

from t 



Chevalier fc Li 



le.g. 



199^ 



200' 



The frequency is inde pendent of s, while Fm oc 1^/^^+"^^) 



^; Sari. Piran. fc Narayan 1998 ). The implied increase in density 
30 s to t 10^ s is a factor of ~5. Such a picture could potentially describe the 
prompt emission via the same external shock that later generates the afterglow. An exter- 
nal shock explanation for the pro mpt emission is motivated superficially by the apparent 



smoothness of the GRB (see, e.g.. iMcMahon. Kumar, fc Panaitescull2004j ). Minimally, the 



model must be able to reproduce the ext remely hard GRB sp e ctrum , with — ^ 0.3 — 0.5 



extending to ~ 1 MeV. During the iBlandford fc McKed (119761 ) evolution stage when 
the fireball expands adiabatically, neither the synchrotron peak frequency {I'm) nor the syn- 
chrotron cooling frequency (z/^) can decrease more rapidly than t"^/^ and go from such a high 
value at t 10 s to the derived values above at t ~ 10^ s. However, prior to this, during the 
fireball deceleration phase when the Lorentz factor F is expected to be roughly constant, the 
cooling frequency will decrease very rapidl y, Ur, oc Our fit s of this model, however. 



with an initial Lorentz factor Fq ~ 700 (e.g.. lMeszaros fc Reeslll997l ) for t < td = 10 s (Figure 



[H]), show that the external shocks which generate the optical fiux underproduces gamma-rays 
by a factor ~ 10. Sync hrotron self-compt on emission cannot help due to that mechanism's 



inefiiciency at small (|Sari fc Esin 



fiux from the IBlandford fc McKed ( 



2001). A more caref ul calculation of the expected GRB 



19761 ) solution (e.g.. ICranot. Piran. fc Sarilll999l ). with 



modifications to the standard synchrotron spectrum to possibly account for the harder GRB 
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versus optical spectrum, is beyond the scope of this paper. 



4.4. Energy Injection 



Refreshed shocks offer a plausible and perhaps better-travelled explanation (e.g. JNousek et al 



20061 ) for the flat early light curve. FigureOshows the external shocks model discussed above, 



but with the addition of a changing shock energy E oc t"- for t < tbreak- For a slow-cooling. 



(Sari & Meszaros 


2000), where a = 


(e.g., 


Sari. Piran. & Naravan 


1998) 



3(1 — p)/4: is the decay index without energy injection 
We derive a = 0.20 ± 0.05. Between 60 and 700 s, the 
0. If the energy iniection is due to a changing luminosity of 



shock energy increases by ^ 

the long-lasting central eng i ne (jllees fc Meszaros 2000; MacFadyen et al. 2001 



Ramirez-Ruiz 



2004J : iLee fc Ramirez- RuizI l2002l ). t he luminosity goes as which may be challeng- 

ing for the progenitor models (e.g., iMacFadyen et al.ll200ll ). If instead the central engine 
generated a flow with a powerlaw distribution of Lorentz factors F, with the slower shells 
of material gradually catching up to the shock, we find M(F) oc r -^-^^*""-^ ( iRees fc Meszaros 
19981 : ISari fc Meszaroslboool : iRamirez-Ruiz. Merloni, fc ReeslboOlh . 



In the constant-density case, the optical flux is more sensitive to the energy injec tion 
than is the X-ray flux (by a factor Aa = a/4; e.g.. Equation (2) in lNousek et al.l 120061 ). If 
X-ray data were available for GRB 051111 during the energy injection episode, we could 
therefore verify and possibly better test the scenario. Long wavelength observations are 
particularly important for testing t he claims of strong energy injection episodes (Aa ^ 1), 
inferred from the X-ray data alone (INousek et al.ll2006l ). 



4.5. The Optical Flux as Reprocessed Prompt Emission 



Similar to what was found above for the increasing density model, the GRB flux at 15 
keV in the energy injection scenarios is ~ 10 times brighter than the expected flux from the 
external shock (solid line in Figure E]). There is a possible steep ening of the light curve decay 
near t f^^: 60 s apparent from ROTSE (IRujopakarn et al.ll2005l ). Super-LOTIS ( iMilne et al. 



20051 ). and the first KAIT observation (Figured]). This early flux decline is too steep, even 
assuming that energy injection has not yet begun to occur at this epoch. These three data 
points therefore suggest that the earliest optical flux is dominated by emission from the 
GRB. 



The dotted curve in Figure [6] shows a powerlaw flt to the gamma-ray data overplotted 
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on the optical. The imphed broadband slope is /3opt-7 = —0.7, consistent with the afterglow 
spectral slope in the optical but marginally inconsistent with the prompt gamma-ray spectral 
slope (Section 13. 2p . If there i s a time delay for the reprocessing of the gamma-ray photons as 
has recently been proposed (IVestrand et al.ll2006l ). this is a coincidence and the broadband 
slope could easily be consistent with the gamma-ray spectral slope. In any case, for this and 
the rising density model discussed above, the synchrotron peak frequency passes through the 
optical quite early {t ^60 s), implying a very inefficient transfer of the shock energy to the 
synchrotron emitting electrons (Table 3), which may also be true for the prompt emission. 



4.6. Light Curve Variability 



The optical light curve at t ;^ 10"^ s (subpanel of Figure [21 Section 13.11) and the opti- 
cal/IR/UV light curve in the range 5 x 10'^ — 1.2 x 10^ s show evidence for residual variability 
at the ^ 30% level, with dt/t ~ 0.1. There is little evidence for spectral change during the 
variability. Such variability may be common and h as been seen previousl y for well sampled 
GRB afterglows. The light curve of GRB 021004 



prominent bumps (IBersier et al.ll2003 



Mirabal et al 



sampled light curve of GRB 030329 (ILipkin et al.l l2004l : IVanderspek et al.l l2003l ) displayed 



Shirasaki et al 



2003 



Fox et al. 



2002h displayed several 



2003a[) . The exquisitely 



prominent departures from a broken powerlaw fa de, even when the underlying supernova 
emission was subtracted away (IBersier et al.ll2003l ). 



Nakar. Piran. fc GranotI (120031 ) explore several possible explanations for the variability 



in GRB 021004 and find that an explanation in terms of in terms of refreshed shocks, a non 



2000; 


Dai & Lu 


2000; 


Lazzati et al. 


2004) 



carry out a similar analysis for GRB 030329 and suggest that density variations are unlikely 
to be responsible for the optical variabili ty due to the appar ent passage of the synchrotron 



cooling break through the passband (e.g., IBersier et al.ll2003l ). Late-time variability after an 



apparent jet break also argues against the patchy jet model, although that model could still 
describe the variability observed at earlier times. 

For GRB 051111, the variability occurs during the phase {t ^ 700 s) when the afterglow 
appears to be well described by external shocks. We believe that this favors an explanation 
in terms of density variation in the surrounding medium. Because the X-ray data are not 
expected to be affected by small variations in the density, this explanation is backed up by 
the quality of the X-r ay fit with a sin g le tem poral powerlaw, although the X-ray error bars 
are ~ 30%. Recently, iGuidorzi et al.l (120051 ) argue that density enhancements can explain 
an achromatic optical light curve bump in data taken 3 minutes after GRB 050502A. More 
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finely sampled data, with broad spectral coverage, will be required to pin down the true 
source of this variability in GRB afterglows. 



Conclusions 



We have presented a thorough analysis of the rich broadband data available for Swift 
GRB 051111, with a focus on the early, multi-band optical data from KAIT and Lulin. 
The optical data prior to t ^ 700 s show a very fiat decline, with little or no evidence for 
a color change as compared to the data after t ^ 700 s. The data at t ;^ 1 minute to 
several hours after the GRB are well fit using a simple, modified external shock model with 
absorption by gas and dust. The modelling entails energy injection or a rapidly increasing 
density profile with radius, prior to the break time. At later times there are large 30%), 
possibly achromatic modulations in the optical, IR, and UV about the best-fit model. The 
modulations appear to be common in well-studied optical afterglows, yet their origin remains 
mysterious. 

The increasing density model may allow for an external shock explanation of both 
the prompt gamma-ray and later emission. Such an increase in density might be ex- 
pected if a supernov a (SN) occurred prior to the GRB, as in the "supra-nova" model 



( Vietri &: Stella IQQsl ). However, the simultaneity of GRB and SN in the nearby, well-studied 
cases (GRB 060218/SN 2006aj, GRB 980425/SN 1998bw, GRB 031203/SN 20031w) argues 
against this possibility. 

On the other hand, slowly declining light c urves in the X-ray band are common and 
are thought to be due to shock refreshment (e.g., iNousek et al.ll2006l ). and the GRB 051111 
observations from KAIT show that this emission can also dominate the optical light curve at 
early times. In fact, the optical light curve for this event, which shows a rapid early decline, 
followed by a levelling off and then a moderate decline typical of those found in the past and 
modelled with external shocks, ap pears quite si r nilar to the "canonical" behavior observed 
in the X-ray band and reported by lNousek et al.l (120061 ). Perhaps there is a canonical optical 
afterglo w behavior, too. (As a counter-point, the early optical behavior reported f or GRBs 
060206 jMonfardini et allbood : Istanek et al]l2006[ l and 060210 Jstanek et al.ll2006h appears 
quite different from that here.) 

An important feature of the GRB 051111 optical afterglow is the lack of a turnover in 
the optical decay rate at early times. We do not detect the peak in the synchrotron spectrum 
passing through the optical bands. Consonant with a lower than average X-ray fiux at late 
times, this leads to a low value for the fraction (assuming constant equipartition) of shock 
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kinetic energy winding up in the synchrotron-emitting elect rons, <; 0.3%. Such a low value 
for ee is uncommon, but not unheard of, in GRBs (see, e.g.. |Panaitescu &: Kumarll2002l ). The 
very low value may indicate a wider diversi ty than previously suspected in the microphysical 
parameters from GRB to GRB (see, also, iBerger. Kulkarni. fc Frailll2003l ). Aside from the 
low eg value, the derived external shocks parameters are comparable to those previously 
found. 

It is important to stress that whether we find ~ 0.3% or ^ 3% depends on the 
deconvolution of the early data into prompt internal shock and afterglow external shock 
components. Our favored smaller eg arises if the synhrotron peak frequency has passed 
through the optical well before the break at t ~ 700 s, in which case mild energy re-injection 
into the shock can explain the gradually decaying light curve. To get the larger eg value, 
the light curve prior to t ~ 700 s must be dominated by the GRB and not by the external 
shocks which dominate after t ~ 700 s. Otherwise, the early light curve would be rising 
rather than declining. In any case, reverbe rations of the prompt emission reprocessed into 
optical radiation (e.g.. lVestrand et al.ll2006l ) appear to be required to explain the earliest few 
optical data points f or GRB 051111. It does not appea r to be necessary to invoke reverse- 
shock emission (e.g.. lAkerlof et al.lll999l : iLi et al.ll2003bl ). which is expected to produce time 
decays more rapid than those observed. 

The absorption at IR, optical, and UV wavelengths in the observer frame is well fit by 
an SMC extinction curve. Little evidence for the 2175 A dust "bump" and an excellent 
fit of the SMC extinction profile are common features in optical GRB afterglow spectr a 



Vreeswiik et al.ll2004J : Ijakobsson et al.ll2003l : ISavadio k Fallll2004J : IWatson et al.ll2006r ) . 



(e.g. 

Combined with the soft X-ray absorption measurement, there is an implied low dust-to- gas 
ratio and a possible overabundance of t he light-metals relat ive to the Fe-group metals. This 
has also been observed for GRB 050401 (jWatson et al.ll2006l ). The light metal overabundance 
works against a direct associa tion with an SMC-like environment, because the SMC has a 
metalicity ~ 1/10 times solar (jPei 1992 ). 



The SED fitting also implies a low ratio of total to selective extinction, Ry f» 2. 
This is a clue that the absorbing medium exhibits unusual dust properties. The work of 
Galama &: Wijerd (1200 ll ) establishing large typical Nu/Ay values has led to suspicions that 



GRBs could destroy dust out to distance s -R ~ 20 pc (see also. 



Waxmann fc Draine 2000 



Fruchter et al. 2001; Draine fc Hao 2002; Perna fc Lazzati 2002 



Perna et all 120031 ). The 



small dust grains are preferentially destroyed, creating a flat (or "gray" ) extinction curve 
(ICalama k WiiersI I2OO1I : ISavadio k Fahl 120041 : Istratta et al1l2005[ l. A small Ry, however, 
suggests small dust grains are dominant, and our SED is very curved. These facts argue 
against the GRB playing a direct role in deflning the extinction properties. The unchanging 
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optical color implies that the absorbing column is not local to the GRB (i.e., R ^ 0.1 pc) 
and may be associated with a nearby giant molecular cloud or the GRB host galaxy. Time- 
resolved spectroscopy in the optical and X-ray bands is of utmost importance for answering 
these questions. 

For GRB 051111, as for a few other Swift events for which ground-based observers have 

been fortunate enough to capture early data, there is an emerging complicated interplay 
between the GRB and the subsequent shocking of the external medium. It is critical that 
more long-wavelength data be taken and published for other GRBs in order to facilitate 
modelling similar to that performed above. Only with such a disentangling of the competing 
emission processes can we hope to answer open questions regarding which shock components 
are truly the most important and what microphysical parameters define the shocks and 
characterize the transition from internal to external shocks. 
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Table 1: KAIT Photometry of GRB 051111 



t [s] 


Exn fsl 


Mag. 




Band 


43.7 


15 


13.60 


0.02 


R 


73.7 


15 


14.54 


0.02 


V 


QQ 7 




1 3 ^7 

L<J .'J t 


n? 

yj .\J ^ 


\ 


124.7 




14.43 


09 

yj .\J ^ 


vt 

-LV. 


156.7 


45 


15.16 


0.03 


V 


212.7 


45 


14.21 


0.02 


I 


268.7 


45 


15.08 


0.02 


R 


324.7 


60 


15.66 


0.02 


V 


395.7 


60 


14.69 


0.02 


I 


470.7 


60 


15.51 


0.02 


R 


578.7 


120 


15.04 


0.02 


I 


734.7 


120 


15.94 


0.02 


R 


892.7 


240 


15.49 


0.02 


I 


1167.7 


240 


16.40 


0.03 


R 


1448.7 


360 


16.05 


0.03 


I 


1843.7 


360 


16.86 


0.03 


R 


2873.7 


120 


17.27 


0.05 


R 


3026.7 


120 


17.40 


0.06 


R 


3178.7 


120 


17.36 


0.05 


R 


3330.7 


120 


17.41 


0.06 


R 


3486.7 


120 


17.49 


0.06 


R 


3639.7 


120 


17.46 


0.07 


R 


3791.7 


120 


17.51 


0.06 


R 


3944.7 


120 


17.62 


0.07 


R 



Notes: R band magnitudes determined for unfiltered observations (Section [2]). We denote with t 
the elapsed time since the BAT trigger. 
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Table 2: Lulin Photometry of GRB 051111 



t [hr] Exp. [s] Mag. (JMag. Band 



4.2201 


300 


19, 


.52 


0, 


.05 


V 


4.3154 


300 


19, 


.12 


0, 


.04 


R 


4.4107 


300 


18, 


.45 


0, 


.04 


I 


4.5071 


300 


20, 


.38 





.07 


B 


4.6021 


300 


19, 


.56 


0, 


.05 


V 


4.6974 


300 


19, 


.15 


0, 


.03 


R 


4.8012 


300 


18, 


.51 


0, 


.04 


I 


4.9154 


300 


20, 


.30 


0, 


.07 


B 


5.0085 


300 


20, 


.37 


0, 


.08 


B 


5.1018 


300 


20, 


.25 


0, 


.07 


B 


5.3204 


300 


20, 


.32 


0, 


.07 


B 


5.4157 


300 


19, 


.54 


0, 


.04 


V 


5.5110 


300 


19, 


.20 


0, 


.03 


R 


6.0479 


300 


19, 


.76 


0, 


.05 


V 


6.1435 


300 


19, 


.27 


0, 


.04 


R 


6.2387 


300 


18, 


.59 


0, 


.04 


I 


6.4037 


300 


20, 


.55 


0, 


.06 


B 


6.4990 


300 


19, 


.89 


0, 


.06 


V 


6.6893 


300 


19, 


.73 


0, 


.05 


V 


6.7846 


300 


20, 


.88 


0, 


.08 


B 


6.8796 


300 


19, 


.86 


0, 


.05 


V 


6.9813 


300 


19, 


.47 


0, 


.05 


R 


7.0762 


300 


18, 


.68 


0, 


.03 


I 


7.1799 


300 


20, 


.86 


0, 


.08 


B 


7 9790 


ouu 


9n 


7"? 


u. 




1 ) 


7.3662 


300 


20, 


.88 


0, 


.10 


B 


7.5546 


300 


19, 


.94 


0, 


.07 


V 


7.6479 


300 


19, 


.96 


0, 


.05 


V 


7.7926 


300 


19, 


.50 


0, 


.03 


R 


7.8879 


300 


18, 


.78 


0, 


.04 


I 


8.0488 


300 


20, 


.72 


0, 


.08 


B 


8.1421 


300 


21, 


.01 


0, 


.14 


B 


8.2365 


300 


21, 


.03 


0, 


,14 


B 


8.3318 


300 


20, 


.25 


0, 


.10 


V 


8.4257 


300 


20, 


,23 


0, 


.08 


V 


8.5193 


300 


20, 


.18 


0, 


.08 


V 


8.6235 


300 


19, 


.82 


0, 


.06 


R 


8.7190 


300 


19, 


.21 


0, 


.06 


I 



Notes: We denote with t the elapsed time since the BAT trigger. 
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Table 3: Broadband Fit Parameters, Const ant- Density ISM 



Paramater 


Value 


E ■ 


6.2 X 10^2 (fixed) 


-C^lum 


3.46 X 10^8 cm (fixed) 


Z 


1.55 (fixed) 


V-y 


(2.1 ±0.3) (t™,opt/[700 s])i-39±o.05 


P 


2.35 ±0.05 


Ce 


(0.03 ±0.01) (t„,opt/[700 s])i-i2±o.05 


ee 


(0.02 ±0.01) (t™,opt/[700 s])-0-ii±0-i5 


n 


(0.8 ±0.5) (t„,opt/[700 s])0-86±0.12 g^-3 



Notes: The fit in Figure H] uses im,,opt — 700 s for tlie passage of tlie synclirotron pealc frequency 
through the optical. 
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16 - 
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22 -_ 



24 



ai =-0.77 ±0.02 



lO^s 



10^ s 



F — r 



20h^F> ^ 
> 



DC _ 



F — r 



KAIT 
Lulin 
GCN 



a2 = - 0.96 ± 0.04 




I. 



10' 



lO'' 10^ 10 

Time Since BAT Trigger [s] 



10^ 



Fig. 1. — R band light curve and {V — R), {R — I) colors from KAIT and Lulin. The best-fit 



from the GCN ( 


Garimella et al. 


2005: 


Milne et al. 


2005: 


Nanni et al. 


2005: 


RuioDakarn et al. 


2005: 


Rvkoff et al. 


2005: 


Sharapov et al. 


2005: 


Smith et al. 2005 


)■ 
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Fig. 2. — Multi-color light curve of GRB 051111 from KAIT and Lulin. The fit residuals for the 
Lulin data are plotted in the inset. 
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Frequency [Hz] 



Fig. 3. — Broadband model overview, data shown at t = 100 s. The afterglow data in the radio, 
optical, IR, UV, and X-rays are well fit by a external shock expanding into a uniform-density 
medium and emitting synchrotron radiation (Sections 14.11 14.31) . The optical and X-ray data drop 
away from the model curve at high and low energies, respectively, due to an appreciable amount of 
absorption by dust and metals (respectively). Also shown are the prompt gamma-ray spectra from 
Swift BAT and Suzaku WAM. 
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1000.00 



100.00 



>^ 10.00 
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KAIT VR 
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Fig. 4. — The data for t ^ 700 s are well fit by a synchrotron external shock model. The 
optical, UV, and IR data are scaled to R band using the best-fit synchrotron shock plus absorp- 
tion model (Table 3, Ar obs = 0.72 mag, Figure [5|). All data include a 10% systematic error 
added in quadrature to the statistical error. The dip near 10^ s is seen in the R band data 



data (Frail et al. 


2005; 


Garimella et al. 


2005; 


Milne et al. 


2005; 


Nanni et al. 


2005 


Poole et al. 


2005; 


Rujopakarn et al 


2005 


RvkofF et al. 


200d: 


Sharapov et al. 2005; 


Smith et al. 


20051. in addition to 



the data from KAIT, PAIRITEL, L ulin, and Swift. All data are corrected for Galactic extinction 
dSchlegel. Finkbeiner. Davislll998l ). 
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Observed Frequency [Hz] ^ 



3 
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SMC 
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um2 



Rest-Frame Frequency [Hz] 



Fig. 5. — (a) Spectral energy distribution for the optical, IR, and UV data interpo- 
lated to T = 5200s after GRB 051111. All data are corrected for Galactic extinction 
(jSchlegel. Finkbeiner. &: David Il998l ). The dashed curve is a powerlaw with (energy) index 



/3 = —0.68, the best-fit value found from the synch rotron sho ck modelling. The other curves 
represent fits of the empirical extinction models from iPeil (|l992l ). with the models labeled in the 
legend and disc ussed in Section l42l (b; inset) The posterior probability for Ay and Ry, using the 
ReichartI (|200ll ) model prior. The la and 2a confidence contours are overplotted. 
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Time Since Burst [seconds] 



Fig. 6. — The early optical and prompt gamma-ray data fit with two variations on the external 
shock model: energy injection prior to t = 700 s with -Bshock t^''^ (solid lines), and an increasing 
density n (x R'^ prior to t = 700 s (dashed lines). Due to an assumed constant initial Lorentz 
factor To = 700, all curves rise prior to the estimated deceleration time = 10 s. Neither model 
reproduces the strong gamma-ray flux. The best-fit powerlaw for the prompt light curve decline 
{F^ oc t-i-4±o.i) is shown as a dotted line, and also scaled by a factor of 350 to match the first 
three optical data points. Here, we assume that the synchrotron peak frequency crosses the R band 
central frequency at tn = 70s (e.g., Table 3), which is an upper limit. We include the GCN data 



(Milne et al 



2005; 



Rujopakarn et al.ll2005l\ in addition to the data from K AIT and Swift. All data 



are corrected for Galactic extinction (ISchlegel. Finkbeiner. &: DavisI Il998l ). The optical data are 
shown mapped in frequency to the R band using the best-fit shock plus absorption model (Table 
3, Ar obs = 0.72 mag, Figure [5]). 



